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Abstract 

Transport of overdamped particles in an asymmetrically periodic potential is investigated in the 
presence of Levy noise and ac-driving forces. The group velocity is used to measure the transport 
driven by the nonthermal noise. It is found that the Levy flights and ac-driving forces are the two 
different driving factors that can break thermodynamical equilibrium. The competition between 
the two factors will induce some peculiar phenomena. For example, multiple transport reversals 
occur on changing the noise intensity. Additionally, we also find that the group velocity as a 
function of Levy index is nonmonotonic for small values of the noise intensity. 
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I. INTRODUCTION 



In systems possessing spatial or dynamical symmetry breaking, Brownian motion com- 
bined with unbiased external input signals, deterministic or random alike, can assist directed 
motion of particles at submicron scales [1]. This subject was motivated by the challenge to 



explain unidirectional transport in biological systems[2], as well as their potential techno- 
logical applications ranging from classical non-equilibrium models[3] to quantum systems^]. 
Ratchets have been proposed to model the unidirectional motion driven by zero-mean non- 
equilibrium fluctuations. Broadly speaking, ratchet devices fall into three categories depend- 
ing on how the app 

flashing ratchets[6|, and correlation ratchets (7|. Additionally, entropic ratchets, in which 
Brownian particles move in a confined structure, instead of apotential, were also extensively 
studied These ratchets demand three key ingredients [s[ which are (a)nonlinearity: it 
is necessary since the system will produce a zero mean output from zero-mean input in a 
linear system; (b) asymmetry (spatial and/or temporal): it can violate the symmetry of the 
response; (c) fluctuating input zero-mean force: it should break t her mo dynamical equilib- 
rium. 

Most studies on ratchets have referred to the consideration of normal diffusion driven by 
Gaussian noises. However, in the past few years, anomalous diffusion has attracted growing 
attention, being observed in various field of physics and related sciences 
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131 ] . Because 



the Levy flights do not possess a finite mean-square displacement, their physical significance 
has been ignored for a long time. However, in the recent years, the growing experimental 
evidences suggest that there is a need to consider a more general type of noise than Gaussian, 
i. e., Levy noise. Description of physical models in terms of Levy flights becomes more and 
more popular 10|-|23|. They are actually observed in various real s yste ms and are used to 



model a variety of processes such as bulk mediated surface diffusion 



17|, exciton and charge 



le systems or 



transport in polymers under conformational motion [18], transport in mice 
heterogeneous rocks [19], two-dimensional rotating flow [20], and many others 

Very few studies on ratchets have focused on the Levy flights. Recently, Dybiec and 
coworkers 21] studied the minimal setup for a Levy ratchet and found that due to the non- 



thermal character of the Levy noise, the net current can be obtained even in the absence 
of whatever additional time-dependent forces. Del-Castillo-Negrete and coworkers [23| also 
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found the similar results in constant force-driven Levy ratchet for 1 < a < 2 (a is Levy 
index). Rosa and Beims [24J studied the optimal transport and its relation to superdiffu- 
sive transport and Levy walks for Brownian Particles in ratchet potential in the presence 
of modulated environment and external oscillating forces. In these studies the Levy noise 
is an intrinsic driving factor to obtain the net transport. In the classical forced thermal 
ratchets jl, 5|, the driving factor is usually the external ac-driving force. However, what's 
the difference between the intrinsic driver and the external one ? How do Levy flights com- 
pete with the ac-driving forces ? In order to answer these questions, in the present paper, 
we studied the transport of overdamped particles in an asymmetrically periodic potential 
in the presence of the Levy flights and ac-driving forces. Our emphasis is on finding the 
difference between the two driving factors and how the competition between them affects 
the transport. 



II. MODEL AND METHODS 

In this study, we consider the transport of Brownian particles moving in an asymmet- 
rically periodic potential in the presence of ac-driving forces and Levy-stable noises. The 
overdamped dynamics can be described by the following Langevin equation in the dimen- 
sionless form 

^- = -U'{x)+A sm{ut) + ( a {t), (1) 

where A and u are the driving amplitude and frequency, respectively. The prime stands 
for differentiation over x. U(x) is an asymmetrically periodic potential (see Fig. 1) 
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FIG. 1: Asymmetrically periodic potential defined by Eq. (2). 
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U(x) = -£/o[sin(x) + j sin(2x)], (2) 

where Uq denotes the height of the potential and A is its asymmetric parameter. 

( a {t) is white, symmetric Levy-stable noise with independent increments distributed ac- 
cording to the stable density with the index a. The time integral of the Levy noise over an 
increment At 

rt+At , , 

L a>a {At) = J (a(t )dt (3) 



is an a-stable process with stationary independent increments and its characteristic function 



(Fourier transform of probability density function) 10Nl2| 



P L {k, At) = exp(-cr|A;| a At), (4) 

where a is the intensity of the Levy noise and k is wave number. The parameter a G (0, 2] 
denotes the stability index, yielding the asymptotic long tail power law for the ^-distribution, 
which is of the |C| _1_a type. For the special case a = 2.0, i.e., for a Gaussian noise, we 
are led back to the Brownian case. From Eqs. (3) and (4) we can obtain the discrete time 
representation of Eq. (1) for sufficiently small time step At 

x{t n+1 ) = x(t n ) - U'(x(t n ))At + A sin(wt n )At + {2D At) «C«(n), (5) 

where n = 0, 1,2... and ( a {n) is a random number possessing Levy stable distribution. In 
order to compare with the classical forced thermal ratchets, we use D = a/2 to describe 
he noise intensity. In numerical simulations, the corresponding generator is taken from 
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sin(aVQ cos([l - a]V) ^ 

Un) = J^vy^ 1 w ] a ' (6) 

where V is random number uniformly distributed on the interval (— tt/2, tt/2), W is an 
independent random variable distributed according to the exponential distribution with 
unit mean. 

In this study, we mainly focus on the transport of the driven particles. In the Gaussian 
noise driving ratchets, due to the existence of the mean value of the noise, the displacement 
of particles also possesses a mean and the transport can be characterized by the average 
velocity. However, for the noise with distribution of a Levy-stable law with < a < 1, 
the mean of the noise and the overall displacement does not exist. The main feature of 
this distribution is that the tails cannot be cut off, or in other words, rare but large events 
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cannot be neglected. As a consequence, the classical stochastic theory (average velocity), 
which is based on the ordinary central limit theorem, is no longer valid. Recently, Dybiec 
and coworkers 21| proposed a different approach to the Levy ratchet problem based on the 
group velocity analysis for < a < 2. Throughout the paper, we will use this method to 
describe the transport of the particles. 

Median line is a very useful tool for investigation of the overall motion of the probability 



density of finding a particle in the vicinity of x [TJJ. A median line for a stochastic process 
x{t) is a function of qo.^{t) given by the relationship Pr(x(t) < 50.5 (£)) = 0.5. Therefore, one 
can use the derivative of the median to define the group velocity of the particle packet [211]. 

v 9 (t) = (7) 

and this definition is valid even for the case of lacking average current. In the following, we 
mainly focus on the study of the long time group velocity, 

V g = hm (8) 

In our simulations, we have considered more than 10 5 realizations to obtain the accurate 
median. In order to provide the requested accuracy of the system dynamics time step was 
chosen to be smaller than 10 -3 . We have checked that these are sufficient for the system to 
obtain consistent results. 



III. NUMERICAL RESULTS AND DISCUSSION 

Our emphasis is on finding the median and group velocity with definitions in Eq. (7). 
In order to investigate the effects of the interplay between the ac-driving forces and Levy 
flights we carried out extensively numerical simulations. For simplicity we set Uq = 1.0 and 
A = 1.0 throughout the work. 

Fig. 2(a) shows the group velocity V g as a function of the driving frequency u. For 
the case of a = 2.0, the Levy ratchet reduces to the classical forced thermal ratchet l|, |5|. 
In this case, the only resource driving particle current across the barrier is the ac-driving 
force. In the adiabatic limit co — > 0, the ac-driving force can be expressed by two opposite 
forces Aq and — Aq. The particles get enough time to cross both side from the minima of 
the potential. It is easier for particles to move toward the slanted side than toward the 
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FIG. 2: (a)Group velocity V g as a function of the driving frequency uj for different values of Levy 
index a at D = 0.1 and Aq = 1.0. (b)Time dependence of the location of the median qo.§{t) for 
different values of driving frequency 00 at D = 0.1 and Aq = 1.0. 

steeper side, so the group velocity is positive. On increasing the frequency u>, due to the 
high frequency, the particles in one period get more enough time to climb the barrier from 
the steeper side than from the slanted side, resulting in negative group velocity. When the 
ac-driving forces oscillate very fast, the particles will experience a time average constant 
force F = j w F(t)dt = 0, so the group velocity goes to zero. At some intermediate values 
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of co, the group velocity crosses zero and subsequently reverse its direction. However, when 
< a < 2, due to the nonthermal character of Levy noise, Levy noise becomes another 
source driving the particle current. As the index a decreases, the positive group velocity 
decreases while the negative group velocity increases. Especially, for a = 0.7, the Levy 
flights dominate the transport and the group velocity is even always negative. It is obvious 
that the transport driven by the Levy flights is opposite to that driven by the ac-driving 
forces for not too small values of the Levy noise intensity. In the same potential, the ac- 
driving forces will induce a positive group velocity, while the Levy flights will give a negative 
group velocity for given index shown in Fig. 2(a). Interestingly, we also found that group 
velocity tends to a negative constant, instead of zero, at fast-driving limit. In this case the 
effects of ac-driving forces disappear and the Levy flights will dominate the transport. In 
Fig. 2(b), we present the time dependence of the location of the median for a given index 
a = 1.5. It is found that the particles exhibit a motion toward the right direction at to — 0.1 
and opposite direction at to = 10.0. 
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FIG. 3: Group velocity V g versus Levy index a with and without ac-driving forces at D = 0.1. 

Figure 3 illustrates the dependence of group velocity V g on the Levy index a with and 
without ac-driving forces. From the figure we can see that the rectified transport can occur 
even the external driving forces are absent. This is due to the nonthermal character of the 
Levy noise which can break the thermaldynamical equilibrium. Interestingly, one can see 
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that the curves demonstrate nonmonotonic behavior for small values of the noise intensity. 
These phenomena can be explained by the interplay between the potential profile and Levy 
flights. Firstly, the particles stay in the minima of the potential awaiting large noise pulses to 
be catapulted out. For not too small values of Levy index, the Levy flights are shorter, and 
the outliers in the Levy noise are smaller, the distance from minima to maxima dominates 
the transport. In this case this distance is shorter from the steeper side (the left side) than 
that from the slanted side (the right side). Consequently, most of the particles are thrown 
out from the steeper side, resulting in negative group velocity. However, for very small 
values of Levy index, the Levy flights are longer, and the outliers in the Levy noise are 
larger. In this case, the slope of the potential dominates the transport. So it is easier for the 
particles moving toward the slanted side than toward the steeper side and the group velocity 
is positive. Therefore, there exists an intermediate value of a at which the group velocity 
takes its extra value and the group velocity as a function of the Levy index is nonmonotonic. 
However, for large values of the noise intensity, the curves demonstrate monotonic behavior. 
This case was reported in Ref.[21]. When a large ac-driving force (A = 1.0) is added, the 
ac-driving force dominates the transport and the group velocity is positive. However, the 
shape of the curve is similar to that without the ac-driving forces. 

In Fig. 4, the group velocity V g is plotted for different values of amplitude A as a 
function of the noise intensity D. When the driving amplitude A Q is small, the Levy flights 
will dominate the transport and the group velocity is negative. On increasing the amplitude 
A Q , the ac-driving forces gradually dominate the transport and the group velocity crosses 
zero and becomes positive, namely, transport reversal occurs. When D — > 0, the group 
velocity tends to zero for small driving fores (A = 0.0, 0.5, 0.65, and 0.75) and goes to a 
finite value for large driving forces (A = 0.85, 0.9, and 1.0). This is due to the fact that 
in the determined ratchets the net current occurs for large driving forces and disappears for 
small driving forces. When D — > oo, the effects of the potential and the ac-driving forces 
disappear, so the group velocity tends to zero. For suitable amplitude A the group velocity 
can change its direction on increasing the noise intensity. Remarkably, multiple transport 
reversals even occur at A = 0.75 (see Fig. 4(b)) and the group velocity reverses its direction 
twice. The intensive competition between the ac-driving forces and the Levy flights leads 
to this peculiar phenomenon. 

We next investigate the role of the amplitude A on the transport. The results are shown 
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FIG. 4: Group velocity V g as a function of the noise intensity D at a = 1.5 and u; = 0.1. (a) 



A =0.0, 0.5, 0.65, 0.85, 0.9, and 1.0. (b) yl =0.70, 0.75, 0.80. 



in Fig. 5. It is found that the group velocity will tend to a negative constant for very small 
driving forces. As we know that in a classical forced thermal ratchet the velocity 

will tend to zero for very small driving forces. The nonthermal character of the Levy noise 
induces this different behavior. For very large values of the amplitude, the influence of the 
potential and the Levy flights will become negligible and the group velocity goes to zero. 
At some intermediate values of amplitude A , the group velocity crosses zero and reverses 



9 



Journal of Statistical Mechanics: Theory and Experiment (2010) P04010 



0.20 



0.15 - 



0.10 



0.05- 



0.00 



-0.05 





-«-D=0.02 




• -O-D=0.10 




-«-D=0.50 




/ \ 

/ o *\ 

f X 


tin nil 

o-o-o-o-o°" 







0.1 



10 



FIG. 5: Dependence of group velocity V g on the driving amplitude Aq for different values of noise 
intensity at a = 1.5 and u = 0.1. 

it direction. 

IV. CONCLUDING REMARKS 



In this paper, we studied the transport of overdamped particles moving in an asymmet- 
rically periodic potential in the presence of ac-driving forces and Levy flights. Because the 
mean of the Levy noise and the mean of the overall displacement do not exist, the definition 
of the average velocity may be invalid. Therefore, throughout the study we use the group 
velocity proposed by Dybiec and coworkers [22] to measure the transport of overdamped 
particles. From the Langevin numerical simulations, we found that ac-driving forces and 
Levy flights are the two different drivers in nature that can break thermodynamical equilib- 
rium. The competition between the two driving factors induces some peculiar phenomena. 
Due to the Levy flights, the group velocity tends to a negative constant, instead of zero, 
for fast-driving limit and it may be always negative for small value of the Levy index. We 
also found that the relation between the group velocity and Levy index is nonmonotonic for 
small values of the noise intensity. There exists an intermediate value of Levy index at which 
the group velocity takes its extreme value. Remarkably, multiple transport reversals occur 
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when the noise intensity changes. This is caused by the competition between the ac-driving 
forces and Levy flights. 

Though the model presented does not pretend to be a realistic model for a real system, 
beyond its intrinsic theoretical interest, the results we have presented have potential appli- 
cations in many processes such as diffusive transport in plasmas, particles separation with 
non-Gaussian diffusion, and ratchet transport in biology systems that are intrinsically out 
of equilibrium. 
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